Next Article in Journal
Multi-Level Alignment Network for Cross-Domain Ship Detection
Previous Article in Journal
A Multi-Objective Quantum Genetic Algorithm for MIMO Radar Waveform Design
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ecological Engineering Projects Shifted the Dominance of Human Activity and Climate Variability on Vegetation Dynamics

1
Key Laboratory of Ecosystem Network Observation and Modeling, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(10), 2386; https://doi.org/10.3390/rs14102386
Submission received: 25 March 2022 / Revised: 13 May 2022 / Accepted: 14 May 2022 / Published: 16 May 2022

Abstract

:
Global greening and its eco-environmental outcomes are getting mounting international focus. The important contribution of China to the global greening is highly appreciated. However, the basic driving forces are still elusive. The Loess Plateau (LP) and Three-River Source Region (TRSR) were chased as study areas in Northern China. The prior one represents the region experiencing intensive human interventions from ecological engineering projects, while the latter is a typical region that is experiencing faster climate change. Hypothesized to be driven by a disproportionate rate of human activities and climates, also being regions of typical large-scale ecological engineering projects, the study goal is to identify the actual driving forces on vegetation dynamics in these two regions. Trend analysis, correlation analysis, and residual trend-based method (RESTREND) were utilized to understand the relationships between climate variability, human activities, and vegetation dynamics. The spatiotemporal variations of vegetation from 1982 to 2019 were evaluated and the respective impacts of climatic and anthropogenic factors on vegetation dynamics were disentangled. Indicating apparent vegetation restoration in LP and TRSR, the results depict that annual LAI has remarkably increased during the 38 years. Temperature and precipitation promoted vegetation growth, whereas the solar radiation and vapor pressure deficit hampered it. After implementing the ecological engineering projects, the primary climatic factor changed from temperature to precipitation. Meanwhile, human activities act as the major driving factor in vegetation greening in the entire study area, with a contribution rate exceeding 70%. This information highlights that ecological engineering can significantly reduce the risks of ecosystem degradation and effectively restore vegetation, especially in ecologically sensitive and vulnerable areas.

Graphical Abstract

1. Introduction

As an indispensable part of the terrestrial ecosystem, vegetation works as a bridge between land and the atmosphere for carbon, water, and energy exchanges, providing food, fiber, fuel, habitat, and other ecosystem services [1,2,3]. Vegetation dynamics are a typical sign measuring plant health and ecosystem stability, which are highly linked to sustainable development [4,5]. Global vegetation is undergoing major changes driven by climate variability and human activities [6,7,8,9]. Climate variability affects vegetation dynamics by affecting their growing environment [10,11,12]. Moreover, human activities such as farming, grazing, and ecological engineering also strongly impact vegetation growth [13,14,15]. Monitoring vegetation dynamics and attributing their underlying driving forces is critical to global vegetation regeneration and conservation in the face of increased human involvement and continued climate change.
Global greening has been plentifully reported. China mainly contributes to global greening, particularly for those regions under the several large-scale ecological engineering projects, such as croplands for grassland and vegetation reforestation on the Loess Plateau (LP) [16]. In the last few decades, the vegetation coverage on the LP has increased by more than 20 percent [17]. A typical combination of human activities and climate variabilities involving vegetation greening within the same type and vegetation type transitions caused this increase. Many studies have concluded that this vegetation greening is mainly caused by human activities [15,18]. Although the impact of ecological engineering on vegetation dynamics seems to be visually clear [18,19], the relative contribution from human activity and climate variability is still exposed to high uncertainties when they are estimated by various methods. To deepen our understanding of the actual driving forces, simultaneously investigating several regions hypothesized to be driven by distinct magnitudes of human activities and climate variabilities is necessary.
Satellite-based remote sensing is an essential means to detect large-scale vegetation dynamics both spatially and temporally compared with ground-based observations [20]. To improve the accuracy of monitoring, selecting favorable remote sensing products is crucial for evaluating vegetation dynamics and their growth conditions. Although the satellite-based normalized difference vegetation index (NDVI) has been largely utilized as a proxy to assess the greenness of the global land surface [21,22], there are several limitations inherent in remote sensing, such as their low sensitivity to vegetation dynamics for dense vegetation cover [23]. Compared with NDVI, the leaf area index (LAI, described as one-half of the total area projected relative to the horizontal size) has a more straightforward physical interpretation. It is also directly connected to the vegetation growth environment [1,24].
Quantitative models have been commonly utilized in recent years to simulate and disentangle the effects of climate variabilities and human activities on vegetation dynamics [25,26]. Frequently used methods include the rain use efficiency (RUE) [27,28], the coefficient of variation method [29], the residual trend (RESTREND) method [30], and the biophysical model-based method [31]. Among them, the RESTREND has been depicted to be an effective approach at large spatial scales [32], and it is commonly accepted due to the practicality, convenience, and universality of its parameters. The RESTREND separates the effects of non-climate factors on vegetation dynamics by calculating actual and simulated vegetation trends based on the relationship model between vegetation and climate factors. However, there are still uncertainties in this method. For example, RESTREND is highly dependent on the established climate-vegetation relationship. Most of the past RESTREND are mainly focused on temperature and precipitation in terms of climates [33] while neglecting some other factors crucial to vegetation growth, such as solar radiation and vapor pressure deficit [34,35,36]. These omissions may significantly impair the effectiveness of the RESTREND.
The LP and Three-River Source Region (TRSR) are fragile and sensitive to climate variabilities and human activities interventions. Significant climate variations have occurred in these areas, which may promote or inhibit vegetation growth [37]. Several studies have revealed that global warming stimulates vegetation growth by extending the growing season and promoting summer photosynthesis. Enhanced precipitation can alleviate drought limitation to some extent, and a higher vapor pressure deficit normally causes inhibitory effects on vegetation growth. Solar radiation can also play a crucial role in vegetation growth [34,35,38,39,40]. On the other hand, on steep slopes overgrazing, deforestation, and other human activities have caused a series of ecological problems over the LP and TRSR, such as vegetation degradation, soil and water loss, and desertification. The Chinese government has invested heavily in launching several large-scale ecological engineering projects to recover and protect environments, such as the Grain for Green (GFG) Project over the LP in 1999 [41] and the Project of Ecological Protection and Construction of the Three-River Source Nature Reserve (EPCP) over the TRSR in 2005 [42]. Propelled by these endeavors, varying degrees of vegetation recovery have occurred in both areas, and land cover in the LP has dramatically changed [43].
Adjacent to the LP, the TRSR is a critical part of the Tibetan Plateau and serves as the water sources region of the three main rivers in China. The ecological engineering projects in the TRSR have been chiefly about grazing for greening, and vegetation coverage change is primarily shaped by natural vegetation growth, which can be hypothesized to be more influenced by natural climates. On the LP, ecological engineering projects consist of a variety of measures, including reforestation and croplands for grassland and forest. Using the TRSR as a reference, we attempt to have a more comprehensive perception of the driving forces on vegetation greening in the two regions.
We will quantify the relative contribution of climate variabilities and human activities to vegetation greening on the LP and TRSR using the RESTREND built on the comprehensive effects of four climate factors (including temperature, precipitation, vapor pressure deficit, and solar radiation). This study aims to: (1) detect the spatial-temporal patterns of vegetation dynamics in the LP and TRSR; (2) explore the relationship between vegetation dynamics and climate variation; (3) quantify the relative contributions of climate variabilities and human activities to vegetation greening and reveal the magnitude differences between LP and TRSR. As global climate variations and human activities intensify, the attribution of vegetation dynamics becomes complex. Improved understanding of the relationships between climate variabilities, human activities, and vegetation dynamics implies utmost importance to climate change mitigation and vegetation restoration.

2. Materials and Methods

2.1. Study Area

Considering the distinct proportion of climate variations and ecological engineering effects on vegetation, we chose two specific regions located in northern-central China as the study area, namely Loess Plateau (LP) and Three-River Source Region (TRSR) (Figure 1). More specifically, LP (100°54′–114°33′E, 33°43′–41°16′N) is situated in the middle reaches of the Yellow River and covers parts of Inner Mongolia, Qinghai, Gansu, Ningxia, Shaanxi, Henan, and Shanxi provinces, with a total area of more than 620,000 km2. The terrain in this region is high in the west and low in the southeast, with numerous ravines and fragmented landforms. The LP is prevailed by an arid and semi-arid climate, transiting from the warm temperate to the cold temperate climate from south to north. The regional climates are characterized as being dry, strong solar radiation and high evaporation in the LP. The TRSR (89°22′–102°26′E, 31°36–37°12′N) is located in the south of Qinghai province with an area of approximately 363,000 km2. This region is the headstream of the Yellow River, the Yangtze River, and the Lantcang River, which is also known as the ‘Chinese water tower’ [42]. The terrain in TRSR rises from east to west, with an altitude range of 2579–6824 m. The TRSR features a typical continental plateau climate characterized as being cold, dry, and intense radiation.

2.2. GLOBMAP LAI Data

The GLOBMAP LAI version 3 dataset (1981–2019) was generated through a fusion of Moderate Resolution Imaging Spectroradiometer (MODIS) LAI and Advanced Very High Resolution Radiometer (AVHRR) LAI data quantitatively, at 8 km resolution and on a Geographic grid [44]. Based on MODIS land surface reflectance data (MOD09A1), the MODIS LAI series was produced by the GLOBCARBON LAI algorithm [45]. The relationships between AVHRR observations (GIMMS NDVI) and MODIS LAI were established pixel by pixel during their overlapped period (2000–2006) [46]. Then the AVHRR LAI back to 1981 was generated by combining historical AVHRR observations and the established pixel-level relationships. Validated by field measurements and fine resolution LAI, this long-term LAI could explain 71% of the variability in the ground LAI over global sites covering all major vegetation types and is biased by 0.81 LAI on average. More detailed descriptions of the algorithm and data evaluation can be found in Liu’s article [44]. Compared with other LAI products, the GLOBMAP LAI exhibits high reliability over China (R2 = 77%) and the low LAI area [47,48]. In addition, this LAI product is the most intra-consistent over time and has a longer time coverage [49]. Its applicability in detecting and attributing terrestrial vegetation greenness has been well-documented at regional and global scales [7,26].
In this study, the GLOBMAP LAI datasets were averaged over the entire year to obtain the annual LAI from 1982 to 2019. We used the annual LAI to reflect vegetation dynamics rather than using growing season LAI, as the prior one is more applicable to our study regions with different growing season lengths. The annual LAI possesses the advantages of being simple and being capable of capturing vegetation greening due to extended vegetation growing season or heightened amplitude of vegetation growth. The annual LAI can further avoid the limitations of growing season LAI, i.e., the subjective thresholds set to define growing season start and end dates [16].

2.3. Meteorological Data

Monthly average temperature (TMP) and total precipitation (PRE) with a spatial resolution of 1 km from 1982 to 2019 were produced by the Loess Plateau science data center, National Earth System Infrastructure of China [50]. Monthly average shortwave down surface radiation (RAD) and vapor pressure deficit (VPD) was obtained from the TerraClimate dataset [51], with a resolution of 4 km. To make them comparable, all meteorological data products were uniformly resampled to the same spatial (8 km) and temporal resolution (annual).

2.4. Auxiliary Data

Digital elevation model (DEM) data with a spatial resolution of 30 m was provided by the Geospatial Data Cloud site, Computer Network Information Center, Chinese Academy of Science (http://www.gscloud.cn (accessed on 28 September 2020)). Land use and cover data in 1980, 2000, 2005, and 2018 were obtained from the Resources and Environment Data Cloud Platform (http://www.resdc.cn (accessed on 30 September 2020)) with a spatial resolution of 1 km. The land was reclassified into cropland, forest, shrubland, grassland, waterbodies, urban and desert.

2.5. Trend Analysis

Trends of LAI, TMP, PRE, RAD, and VPD over the past 38 years were assessed based on simple linear regressions. Ecological engineering projects were launched in 2000 over the LP and 2005 over the TRSR. We split the entire study period (1982–2019) into the first period (1982–1999 for the LP and 1982–2004 for the TRSR) and the second period (2000–2019 for the LP and 2005–2019 for the TRSR) following Xiu [41] and Zhai [52]. The trend slope in a multi-year regression equation represents interannual changes and can be solved using the ordinary least squares method (OLS). The equation is as follows:
β = n × i = 1 n i × x i i = 1 n i i = 1 n x i n × i = 1 n i 2 ( i = 1 n i ) 2
where β represents the interannual trend of x i for each pixel, n is the number of years, and x i represents the value of x in the i th year. β > 0 suggests an increasing trend, while β < 0 suggests a decreasing trend. The significances of the trends were determined using a t-test at a significance level of 0.05.

2.6. Correlation Analysis

Relationships between LAI and climate factors were assessed utilizing simple parametric correlations of the Pearson product-moment linear correlation coefficients ( r values ranging from −1 to 1) as follows:
r = n i = 1 n x i y i i = 1 n x i × i = 1 n y i n i = 1 n x i 2 ( i = 1 n x i ) 2 × n i = 1 n y i 2 ( i = 1 n y i ) 2
where r represents the Pearson correlation coefficient, n is the number of years, x i refers to LAI in the i th year, and y i represents climate factor in the i th year. The significance of the correlations was assessed by t-test at significance levels of 0.05.

2.7. Residual Trend Analysis

Due to the lack of long-term spatially explicit data about human activities such as grazing, mining, and land use management, it is normally assumed that effects of human activities are the residual after climate variabilities effects are accounted for, which can be disentangled by RESTREND [53,54,55]. This study thoroughly considered the complexity of the multi-layer background features. We established the relationship model between vegetation dynamics and climate factors by taking TMP, PRE, RAD, and VPD as independent variables and LAI as the dependent variable. The RESTREND method works as follows:
Step 1: Establish a multiple linear regression between LAI and climate factors as follows:
L A I p r e = A × T + B × P + C × R + D × V + E
where, L A I p r e represents the simulation value of LAI. T represents the TMP, P represents the PRE, R represents the RAD, V represents the VPD. A, B, C, and D are the slopes of the linear regression between the LAI and TMP, PRE, RAD, and VPD, respectively. E is a constant.
Step 2: Calculate the residual difference between L A I a c t u a l and L A I p r e , and define the residual differences as L A I r e s as follows:
L A I r e s = L A I a c t u a l L A I p r e
where L A I r e s represents the residual, L A I a c t u a l represents the actual value of LAI proxied by satellite products, L A I p r e represents the simulation value of LAI.
Step 3: In this study, we calculated the trends of L A I a c t u a l , L A I p r e and L A I r e s using the simple linear regression according to Equation (1). The trends of L A I p r e represent the impact of climate variabilities, and the trends of L A I r e s measure the effects of human activities.
Step 4: Set the division rule of driving factors on vegetation greening as shown in Table 1. Based on this, the contribution of climate variabilities and human activities to vegetation greening before and after the implementation of ecological engineering was calculated.

3. Results

3.1. Spatial-Temporal Patterns of Vegetation Dynamics

The annual average LAI trends from 1982–2019 revealed that vegetation growths in the LP and the TRSR have experienced significant amelioration (Figure 2). Various ecological engineering projects have different implementation timelines., we set the separation line as 1999–2000 for the LP and 2004–2005 for the TRSR. In the LP, the annual LAI depicted a substantial upward trend at a rate of 0.0050 year−1 (p < 0.01) from 1982 to 2019. The LAI trend in 1982–1999 exhibited a non-significant increase (0.0014 year−1, p > 0.05), but it sharply increased in 2000–2019 with a rate of 0.0133 year−1 (p < 0.01). In the TRSR, the background annual LAI is low and showed a fluctuating upward trend during 1982–2019 with a rate of 0.0012 year−1 (p < 0.01). The annual LAI fluctuation in 2005–2019 was more noticeable than in 1982–2004.
Vegetation dynamics in the study areas exhibited high spatial heterogeneity (Figure 3a). The annual LAI in 97.5% of the LP and 86.0% of the TRSR showed increasing trends, with 91.4% and 65.3% of each region showing a significant trend (p < 0.05), respectively (Table 2). The noticeable upward trends (more than 0.01 year−1) occurred mainly in the southeastern LP region dominated by forests. The decreasing trend accounted for 2.5% of the LP and 14.1% of the TRSR. However, the area exhibiting significant decreasing trends was small (only 0.5% of the LP and 3.5% of the TRSR), dispersed in the southern region of LP and the central and south parts of TRSR.
Vegetation trends varied considerably before and after implementing the ecological engineering project (Figure 3b,c). LAI showed an increasing trend in 80.3% of the LP and 66.2% of the TRSR before the implementation, with 23.3% and 26.4% of the total being significantly increased, respectively. The area with major vegetation improvements was mostly distributed in the northern region of LP and the northwestern region of TRSR; 0.9% of the LP and 3.9% of the TRSR suffered a certain degree of vegetation degradation, mainly around the area with rapid urbanization. LAI exhibited a significant increase for almost the entire LP after the implementation, and 51.1% of the LP demonstrated a steeper increasing slope (more than 0.01 year−1). On the other hand, only 17.0% of the TRSR showed a significant increasing trend, and decreasing trends were mainly found in the central and southeastern regions, accounting for 2.6% of the TRSR. During the second period, in the LP, the area of significant increasing trends expanded by 71.4%, and the area exhibiting significant decreasing trends shrank by 3.6%. However, in the TRSR, the area showing significant increasing and decreasing LAI trends shrank by 9.4% and 1.3% after 2004, respectively. Overall, vegetation showed a greening trend during 1982–2019, and the area exhibiting a greening trend expanded after the implementation of the ecological engineering project, especially in the LP.

3.2. Climate Variability in the Study Area

The annual mean temperature across the LP and TRSR has been increasing at a significant pace of 0.0363 °C/year (p < 0.01) and 0.0256 °C/year (p < 0.01) over the past 38 years (Figure 4). The annual total precipitation increased significantly in the LP and exhibited a weak increase in the TRSR. The annual RAD and VPD showed opposite trends between the LP and TRSR during 1982–2019. In the LP, the RAD exhibited non-significant trends, and VPD increased at a rate of 0.0039 kPa/year (p < 0.01). However, in the TRSR, the RAD showed a significantly declining trend with a rate of −1.672 W/m2/year (p < 0.01), and the VPD decreased non-significantly.
The climate variability in the study area exhibited high spatial and temporal heterogeneities (Figure 5). Before implementing the ecological engineering project, almost the whole area showed increasing TMP and VPD in the LP. In contrast, 68.3% of the area showed a decreasing PRE and 41.3% of the area exhibited an increase in RAD. Approximately 70% of the LP showed warming and drying climates during 1982–1999. However, nearly 80% of the LP slowed its warming and became wetting from 2000 to 2019. In the TRSR, the west part showed warming and wetting climates while the east showed warming and drying during 1982–2004. The climate of the whole TRSR then reversed as it transitioned from the first to the second period, as seen by cooling and drying conditions in about 80% of the area.

3.3. The Relationship between Vegetation Dynamics and Climate Variability

Due to the tremendous spatial heterogeneity of the climate variability in the study area, the relationship between vegetation dynamics and climate variability also exhibited substantial spatial and temporal patterns (Figure 6). During 1982–2019, the correlation coefficient between LAI and TMP was positive in 96.65% of the LP and 90.05% of the TRSR, with 72.09% of the LP and 57.09% of the TRSR being significant, respectively (Table 3). The rising TMP had a positive effect on vegetation growth in the whole period. In almost 85% of the study area, after implementing the ecological engineering project, the correlation between LAI and TMP remained positive. The correlation coefficients between PRE and LAI in the majority of the study area (92.74%) were positive during 1982–2019. Compared with the first period, the correlation shifted from negative to positive in the second period. The positive influence of PRE on LAI became stronger. Averaged over the past 38 years, the correlation between LAI and RAD was gradually changed from negative in the west to positive in the east. The RAD had a positive effect on LAI in about 50% of the LP during the three periods, while in the TRSR, the RAD mainly had a negative impact (about 80% of the area) on LAI. Between 1982–2019, the VPD was positively correlated with LAI in the LP (80.24%) and the linkage became negative in the TRSR (73.29%). However, in 83.34% of the LP, the correlation between LAI and VPD changed to negative after implementing the ecological engineering project. Compared with 1982–2004, the positive correlation area of the TRSR between LAI and VPD expanded by 23.6% during 2005–2019.
According to the coefficients of correlation and their significances between LAI and climatic drivers, we draw maps showing the spatial separation of dominant climatic constrain on LAI (Figure 6e1–e3). From 1982 to 2019, temperature was the dominant climatic factor in over 50% of the entire study area. More than 10% of the area is dominated by precipitation, mainly in the eastern part of the TRSR and the western part of the LP. After implementing the ecological engineering project, vegetation growth in the study area changed from temperature to precipitation constraints. The area of RAD dominance increased in the TRSR, and the impact of VPD in the LP was enhanced and spatially more concentrated.

3.4. Land Cover Changes before and after the Implementation of Ecological Engineering

Before implementing ecological engineering, the areas of cropland, forest, and urban on the LP increased, while the areas of other land cover types reduced (Table 4). Grassland was the main source of increase in the area of other land cover types, with 3364.03 km2 of grassland converted to cropland and 1568.61 km2 converted to desert. In the TRSR, the areas of cropland, forest, urban, and desert increased, whereas others decreased. Land cover change in the TRSR was relatively weak before the implementation of ecological engineering, with mainly grassland (924.35 km2) and waterbodies (485.08 km2) converting to desert.
The areas of cropland and desert on the LP significantly decreased after implementing ecological engineering, covering the areas of 11,550 km2 and 4298 km2, respectively. The reduced cropland area was mainly changed into grassland, urban, and forest, whereas the reduced desert area was mainly converted into cropland and grassland. In the TRSR, the areas of grassland and waterbodies increased significantly and the area of desert decreased significantly, accompanied by the non-significant area changes of other land cover types. The area of grassland mainly increased, stemming mainly from the desert. During this period, the waterbodies area increased, stemmed from grassland and desert with a total of 7000 km2. Despite the area of flow between grassland, forest, and shrubland being large, the transfer in and out was basically balanced, so the areas of forest and shrubland remained stable.

3.5. Contributions of Climate and Human Factors to Vegetation Greening

In this study, we quantitatively identified the contributions of climate variabilities and human activities to vegetation greening by using a linear regression model. Figure 7 showed the linear regression coefficients between LAI and TMP, PRE, RAD, and VPD. During 1982–2019, climate variabilities and human activities jointly contribute to vegetation greening (Figure 8). The average contributions of climate variabilities to vegetation greening were 59.91% in the LP and 52.92% in the TRSR, both slightly higher than the contributions of human activities. Moreover, climate variabilities dominated vegetation greening accounting for 81.23% and 52.33% of the LP and TRSR, respectively. On a long-timescale, climate variability contributed more to vegetation greening than human activity.
Nevertheless, the relative contribution of climate variabilities and human activities shifted considerably between the two time periods, divided by the initiation of the ecological engineering projects. Before implementing the projects, climate variabilities controlled 81.52% of vegetation greening area on the LP, and the positive effect of human activities on vegetation was small. In the TRSR, the western side was mainly influenced by climate change and the eastern side by human activities. Furthermore, the impact of climate variabilities (50.58%) was weakly higher than that of human activities. These indicated that vegetation greening in the study regions was mainly dominated by climate variabilities. After the projects’ implementation, the area of vegetation greening mainly caused by human activities increased significantly from 6.69% to 74.62% in the LP and from 27.62% to 93.67% in the TRSR. The patches dominated by human activities (approximately 430,208 km2) were concentrated and large in the LP, while those in the TRSR (approximately 39,744 km2) were small and scattered.

4. Discussion

4.1. Vegetation Greening

This study analyzed spatial-temporal patterns of vegetation dynamics across different periods over the past 38 years, and the characteristics of the LP and TRSR were compared. We detect significant vegetation greening in the LP (91.4%) and TRSR (65.3%) during 1982–2019, with considerable greening mainly distributed in the eastern TRSR and southeastern LP. These findings are in line with earlier studies which show that vegetation conditions of the study regions have improved [15,17,56]. Although LAI decreased in some parts of the study region, the decreasing trend was mostly not statistically significant. The areas where LAI decreased significantly were consistent with the distribution of urban land, indicating that the significant decrease in LAI was mainly related to human activities such as urbanization. However, the distribution of non-significant decreases changed from sporadic to clustered after 2005, reflecting heightened fluctuation of the LAI trend [43]. This suggested that the vegetation of TRSR has generally improved, but vegetation degradation still highly exists.
Although both LP and TRSR exhibit significant greening from 1982–2019, their specific change characteristics differ greatly. First, the rate of vegetation greening on the LP is about four times higher than that of the TRSR, which may be related to therein high background vegetation coverage and the extensive ecological engineering project implementation [42]. According to previous studies, more than 20,000 km2 of cropland has been converted to grassland and forests after the GFG implementation [57,58], as well as the afforestation area on the LP, reached 28,300 km2 in 2018 alone [59]. Second, vegetation dynamics show different trends between the two regions prior to and posterior to the ecological engineering projects. After implementing the ecological engineering projects, the vegetation greening accelerated, and the greening area extended obviously in the LP. By contrast, only 17% of the TRSR exhibited significant greening trends, much lower than 94.7% of the LP. The faster vegetation greening rate on the LP is highly related to the more favorable climate conditions and the types of ecological engineering measures on the LP. On the LP, the combinations of precipitation and temperature produce environmental conditions breeding higher vegetation productivity than on the TRSR. Moreover, ecological engineering projects on the LP are mainly about afforestation and grassland recovery from croplands, which accommodate higher productivity ecosystems than the alpine grasslands on the TRSR.

4.2. Climate Variability and Its Relationship with Vegetation Dynamics

There is no doubt that climate variability is a factor essential for vegetation dynamics [19]. Besides temperature and precipitation, solar radiation and vapor pressure are increasingly included in climate variability and anthropogenic activity differentiation studies [60,61,62]. Our results show that vegetation growth on the LP and TRSR was mainly influenced by temperature, followed by precipitation, vapor pressure deficit, and lastly by solar radiation. For nearly 10% of our study area, vegetation dynamics are primarily influenced by solar radiation and vapor pressure deficit. Our model has significantly improved the variations explained by adding solar radiation and VPD to the statistical model (R2 = 85.3%), as compared to the previous TP-based LAI model (R2 = 65.6%), which considered only temperature and precipitation, our model has significantly improved the variations explained by adding solar radiation and VPD to the statistical model (R2 = 85.3%). Other climate parameters were not included in this analysis, considering the risk of multicollinearity and marginal improvements on explained variations [63,64].
Our study confirms previous reports that vegetation growth is mainly limited by regional hydrothermal conditions in the LP and TRSR [6,36,65]. Temperature rises may promote photosynthesis and extend vegetation growing season, thereby benefiting dry matter accumulation [39,66]. Precipitation favors vegetation growth in both LP and TRSR, which are both classified as semi-regions and enhanced precipitation can alleviate drought limitation to some extent [18,67,68]. Nevertheless, in areas such as southeast LP and central TRSR, precipitation is negatively correlated with vegetation growth. This downplaying effect may be related to the artificial irrigation that makes vegetation not sensitive to precipitation in some parts of the LP [69]. In cold and the relatively humid TRSR, precipitation may decrease temperature and solar radiation, resulting in inhibiting vegetation growth [70]. Temperature and precipitation have similar impacts on vegetation dynamics in most of the studied locations, while the effects of solar radiation and vapor pressure deficit are regionally diverse. On the eastern side of the LP, solar radiation benefits vegetation growth via its promotion of photosynthesis [18,71]. However, in the TRSR and western part of the LP, high-altitude areas with sufficient solar radiation, increased solar radiation leads to more evapotranspiration, thereby inhibiting vegetation growth [72]. The higher VPD normally shows inhibitory effects on vegetation growth because leaf and canopy photosynthetic rates decrease due to stomatal closure [35]. However, VPD and vegetation greening are positively correlated in the LP. The unusual relationship might be caused by the particular soil moisture conditions on the LP. The intervention of human management might be another contributing factor.
The relationship between climates and vegetation adjusted after the ecological engineering projects. In the LP, climate variations shifted from warming-drying to warming-humid after 2000. Although the warming rate becomes slower, precipitation and solar radiation shift from decreasing to increasing in the LP, and vapor pressure deficit changes from increasing to decreasing [41,73]. The LP’s environment becomes more favorable to vegetation growth. However, in the TRSR, warming slowed down after 2005 and the trend was even reversed in some regions. Reduced precipitation, stronger solar radiation, and increased vapor pressure deficit make the climate conditions more unfavorable for vegetation growth in the TRSR [74]. Precipitation has substituted temperature as the dominant climate factor affecting vegetation growth in many regions [75].
In a word, climate variability and its relationship with vegetation dynamics exhibit similarities and disparities between LP and TRSR. After the ecological engineering projects, the dominant climate factor changes from temperature to precipitation in both LP and TRSR, which indicates that human activities might significantly modify the regional hydrological conditions. After the ecological engineering projects, climate conditions in the LP improved while degraded in the TRSR. The progressively unfavorable climate conditions in the TRSR make it more challenging to implement ecological engineering projects in the TRSR.

4.3. Ecological Engineering Drives Vegetation Greening

This study quantifies the relative contribution of climate variabilities and human activities to vegetation greening (p < 0.05) using the RESTREND method. However, both human activities and climate factors are contained in the residual. Some of these climate factors (such as soil texture) are relatively stable on the time dimension, and their impact on vegetation dynamics can be neglected. Human variables are usually assumed to account for a considerable amount of the residual. Furthermore, in the LP and TRSR, most human activities positively impacting vegetation growth are stemmed from ecological engineering projects. Our results identified a significant increase in the contribution of human activities after implementing ecological engineering projects, which are in line with the series of government launched large-scale projects in these two regions. In addition, the changes in land cover before and after the ecological engineering was considered, and the results also confirmed the increased contribution of human activities to vegetation greening.
During 1982–2019, vegetation greening on the LP and TRSR was driven by a combination of climate variabilities and human activities, with the contribution of climate variabilities slightly higher than human activities. This is consistent with prior research findings. [18,52]. However, in the TRSR, the contribution of climate variabilities is lower than human activities in both the two split periods. The explanation for this could be that the relationships between vegetation dynamics and climate are more stable in an extended period, while human activities are exposed to wider fluxes [30,61].
The human activity contribution has been intensified with time elapsing. The GFG project, also known as the Conversion of Cropland to Forests and Grassland Program, was launched in 1999 on the Loess Plateau and has played a significant role in vegetation greening directly or indirectly. The primary measures include converting agricultural lands into forest or grassland, afforestation, artificial irrigation, and so on. From 2000–2014, there were afforested areas of 2.39 × 104 km2 in the northern Shaanxi Province alone. It was reported that the correlation coefficients between the NDVI time series and cumulative afforestation areas in some areas of the LP reached 0.978 [76]. These findings suggest that human activities have been the primary factor in vegetation improvement in the LP. Our results also show that after the GFG project, the greening area primarily driven by human activities on the LP increased from 6.7% to 74.7%. This result reveals that the implementation of ecological restoration projects and the economic incentives given to the local farmers have made a remarkable contribution to vegetation greening.
Intense human intervention has altered the land surface state in the LP, regional climates, and their relationships with vegetation dynamics. Some studies have noticed that the conversions of unused land and rainfed cropland into forest or grassland have contributed to the increasing precipitation trends after 1999 [17]. At the same time, it has raised concerns about local food shortages. However, in Shaanxi Province, the fastest-greening province, the crop production area decreased by 21.34% while the total yield and yield per unit area increased by 12.60% and 43.14% from 2000–2018 [15]. The increased yield per unit area and the optimization of crop structure compensated for the yield loss caused by the reduction of cropland area.
In the TRSR, the first phase (2005–2012) of EPCP comprises reducing livestock, banning grazing, delaying grazing, ecological migration, controlling desertification, rodent pest control, and returning farmland to forests and grasslands [77,78]. The second phase was launched in 2013, and it adheres to the primary principle of reducing human interventions and investing more in natural vegetation recovery. Animal husbandry is the main industry in the TRSR, and the number of livestock has an important impact on vegetation dynamics. Therefore, the change in livestock numbers can reflect the effect of EPCP to some extent. According to the statistical annals of the Qinghai government, livestock numbers in the TRSR have decreased by 15.34% since the implementation of the EPCP [74]. Our results also show that human activities dominated proportion increased from 27.6% to 93.7% for those greening areas on the TRSR. Whereas, climate conditions in the TRSR become increasingly unfavorable for vegetation growth as mentioned above. After the EPCP, vegetation greening is not significant in many areas of the TRSR under the interacted human activities and natural climates. Notably, some protective measures are not as effective as expected. Thus, more appropriate adaptation strategies should be adopted to address climate change in the future.
A high disparity was identified between the LP and TRSR in terms of the contribution rates of climate variabilities and human activities to vegetation greening. Although the vegetation greening rate in the TRSR is smaller than in LP, the ecological engineering project’s contribution rate is higher in the TRSR. Our results, together with other relevant studies, raise the importance of paying more attention to the contribution of artificial vegetation construction and agricultural technology innovation to regional vegetation greening [79,80,81,82]. Thereby, for ecologically fragile areas, vegetation recovery can incline more on natural recovery, but a balanced natural recovery and human interventions are also necessary.
The entire investigation made several improvements and bridged some gaps based on existing studies. Firstly, RESTREND can capture the enhancement of the anthropogenic contribution to vegetation greening after implementing ecological engineering projects in the study areas, and the inclusion of suitable climate factors can make this method more applicable locally. Secondly, by comparing the relationship between vegetation, climate, and human activities before and after the ecological engineering projects, it was unexpectedly found that climate variabilities in the TRSR after 2005 mainly inhibited vegetation growth, and vegetation greening was dominated by human activities. Thirdly, it was easy to overestimate the contribution of human activities because climate effects on vegetation are difficult to detect over shorter time spans. Thus, two typical regions need to be compared in order to have a more accurate assessment of the anthropogenic contribution. Finally, through comparison, it was found that the effects of implementing ecological engineering projects varied greatly in different areas, with the LP exhibiting obvious results while the TRSR exhibited fewer ones. In implementing ecological engineering projects, close attention should be paid to climate variations and timely adjustments to the specific measures of ecological engineering projects, such as increasing human intervention when climate conditions are not suitable for the natural recovery of vegetation.

5. Conclusions

The present study analyzed the spatial-temporal patterns of vegetation dynamics over the LP and TRSR from 1982 to 2019. The relative contributions of climate variabilities and human activities on vegetation greening were quantified using RESTREND based on the LAI-climate model. The following conclusions were drawn: (1) On an overall basis, there was a significant vegetation greening trend in the entire study area. Vegetation greening rate and magnitude were greater on the LP than on the TRSR; (2) Increasing TMP and PRE mainly promote vegetation growth in the entire study area, while higher RAD and VPD mainly play an inhibiting role; (3) Ecological engineering projects have adjusted the relationship between climate and vegetation. The primary climatic factor changed from temperature to precipitation after implementing the series of ecological engineering projects; (4) The most notable vegetation greening primarily stems from the ecological engineering projects. Counterintuitively, the area percentage of human activities dominated region is higher in the TRSR (91.73%) than in the LP (68.02%); and (5) Compared with the Loess Plateau, the results of ecological engineering are more challenging to maintain in the TRSR. This study can further enhance our understanding of the relationship between climate variabilities, human activities, and vegetation dynamics for the large-scale ecological engineering project affected regions.

Author Contributions

Conceptualization, Y.Z. (Yangjian Zhang) and J.G.; Methodology, J.G., Z.Z. and N.C.; Data curation, J.G., Y.Z. (Yixuan Zhu), Y.S. and Y.Z. (Yu Zhang); Formal analysis, J.G., Y.Z. (Yixuan Zhu), Y.C. and G.Z.; Investigation, J.G. and J.Z.; Software, J.G.; Visualization, J.G. and Y.Z. (Yixuan Zhu); Writing—original draft, J.G.; Writing—review and editing, Y.Z. (Yangjian Zhang) and J.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was jointly supported by the National Natural Science Foundation of China (Grant No. 41725003), the National Key Research and Development Program of China (2019YFA0607302), and the National Natural Science Foundation of China (Grant No. 41991234).

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Piao, S.; Wang, X.; Park, T.; Chen, C.; Lian, X.; He, Y.; Bjerke, J.W.; Chen, A.; Ciais, P.; Tømmervik, H.; et al. Characteristics, drivers and feedbacks of global greening. Nat. Rev. Earth. Environ. 2019, 1, 14–27. [Google Scholar] [CrossRef]
  2. Bonan, G.B.; Pollard, D.; Thompson, S.L. Effects of boreal forest vegetation on global climate. Nature 1992, 359, 716–718. [Google Scholar] [CrossRef]
  3. Haberl, H.; Erb, K.H.; Krausmann, F.; Gaube, V.; Bondeau, A.; Plutzar, C.; Gingrich, S.; Lucht, W.; Fischer-Kowalski, M. Quantifying and mapping the human appropriation of net primary production in earth’s terrestrial ecosystems. Proc. Natl. Acad. Sci. USA 2007, 104, 12942–12945. [Google Scholar] [CrossRef] [Green Version]
  4. Yan, W.; Wang, H.; Jiang, C.; Jin, S.; Ai, J.; Sun, O.J. Satellite view of vegetation dynamics and drivers over southwestern China. Ecol. Indic. 2021, 130, 108074. [Google Scholar] [CrossRef]
  5. Li, Y.; Piao, S.L.; Li, L.Z.X.; Chen, A.P.; Wang, X.H.; Ciais, P.; Huang, L.; Lian, X.; Peng, S.S.; Zeng, Z.Z.; et al. Divergent hydrological response to large-scale afforestation and vegetation greening in China. Sci. Adv. 2018, 4, eaar4182. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Ge, W.; Deng, L.; Wang, F.; Han, J. Quantifying the contributions of human activities and climate change to vegetation net primary productivity dynamics in China from 2001 to 2016. Sci. Total Environ. 2021, 773, 145648. [Google Scholar] [CrossRef]
  7. Piao, S.; Yin, G.; Tan, J.; Cheng, L.; Huang, M.; Li, Y.; Liu, R.; Mao, J.; Myneni, R.B.; Peng, S.; et al. Detection and attribution of vegetation greening trend in China over the last 30 years. Glob. Chang. Biol. 2015, 21, 1601–1609. [Google Scholar] [CrossRef]
  8. Sun, W.; Song, X.; Mu, X.; Gao, P.; Wang, F.; Zhao, G. Spatiotemporal vegetation cover variations associated with climate change and ecological restoration in the Loess Plateau. Agric. For. Meteorol. 2015, 209–210, 87–99. [Google Scholar] [CrossRef]
  9. Wang, L.; Tian, F.; Wang, Y.; Wu, Z.; Schurgers, G.; Fensholt, R. Acceleration of global vegetation greenup from combined effects of climate change and human land management. Glob. Chang. Biol. 2018, 24, 5484–5499. [Google Scholar] [CrossRef]
  10. Tao, J.; Xu, T.; Dong, J.; Yu, X.; Jiang, Y.; Zhang, Y.; Huang, K.; Zhu, J.; Dong, J.; Xu, Y.; et al. Elevation-dependent effects of climate change on vegetation greenness in the high mountains of southwest China during 1982–2013. Int. J. Climatol. 2018, 38, 2029–2038. [Google Scholar] [CrossRef]
  11. Chen, B.; Zhang, X.; Tao, J.; Wu, J.; Wang, J.; Shi, P.; Zhang, Y.; Yu, C. The impact of climate change and anthropogenic activities on alpine grassland over the Qinghai-Tibet Plateau. Agric. For. Meteorol. 2014, 189–190, 11–18. [Google Scholar] [CrossRef]
  12. Wang, X.; Wang, T.; Liu, D.; Guo, H.; Huang, H.; Zhao, Y. Moisture-induced greening of the South Asia over the past three decades. Glob. Chang. Biol. 2017, 23, 4995–5005. [Google Scholar] [CrossRef] [PubMed]
  13. Sun, Z.; Mao, Z.; Yang, L.; Liu, Z.; Han, J.; Wanag, H.; He, W. Impacts of climate change and afforestation on vegetation dynamic in the Mu Us Desert, China. Ecol. Indic. 2021, 129, 108020. [Google Scholar] [CrossRef]
  14. Huang, K.; Zhang, Y.; Zhu, J.; Liu, Y.; Zu, J.; Zhang, J. The Influences of Climate Change and Human Activities on Vegetation Dynamics in the Qinghai-Tibet Plateau. Remote Sens. 2016, 8, 876. [Google Scholar] [CrossRef] [Green Version]
  15. Kou, P.; Xu, Q.; Jin, Z.; Yunus, A.P.; Luo, X.; Liu, M. Complex anthropogenic interaction on vegetation greening in the Chinese Loess Plateau. Sci. Total Environ. 2021, 778, 146065. [Google Scholar] [CrossRef] [PubMed]
  16. Chen, C.; Park, T.; Wang, X.; Piao, S.; Xu, B.; Chaturvedi, R.K.; Fuchs, R.; Brovkin, V.; Ciais, P.; Fensholt, R.; et al. China and India lead in greening of the world through land-use management. Nat. Sustain. 2019, 2, 122–129. [Google Scholar] [CrossRef]
  17. Naeem, S.; Zhang, Y.; Zhang, X.; Tian, J.; Abbas, S.; Luo, L.; Meresa, H.K. Both climate and socioeconomic drivers contribute to vegetation greening of the Loess Plateau. Sci. Bull. 2021, 66, 1160–1163. [Google Scholar] [CrossRef]
  18. Li, J.; Peng, S.; Li, Z. Detecting and attributing vegetation changes on China’s Loess Plateau. Agric. For. Meteorol. 2017, 247, 260–270. [Google Scholar] [CrossRef]
  19. Shi, S.; Yu, J.; Wang, F.; Wang, P.; Zhang, Y.; Jin, K. Quantitative contributions of climate change and human activities to vegetation changes over multiple time scales on the Loess Plateau. Sci. Total Environ. 2021, 755, 142419. [Google Scholar] [CrossRef]
  20. Qu, L.; Huang, Y.; Yang, L.; Li, Y. Vegetation Restoration in Response to Climatic and Anthropogenic Changes in the Loess Plateau, China. Chin. Geogr. Sci. 2020, 30, 89–100. [Google Scholar] [CrossRef] [Green Version]
  21. Zhang, Y.; Fan, C.; Huang, K.; Liu, Y.; Zu, J.; Zhu, J. Opportunities and challenges in remote sensing applications to ecosystem ecology. Chin. J. Ecol. 2017, 36, 809–823. [Google Scholar]
  22. Chu, H.; Venevsky, S.; Wu, C.; Wang, M. NDVI-based vegetation dynamics and its response to climate changes at Amur-Heilongjiang River Basin from 1982 to 2015. Sci. Total Environ. 2019, 650, 2051–2062. [Google Scholar] [CrossRef]
  23. Qiao, Y.; Jiang, Y.; Zhang, C. Contribution of karst ecological restoration engineering to vegetation greening in southwest China during recent decade. Ecol. Indic. 2021, 121, 107081. [Google Scholar] [CrossRef]
  24. Lin, J.Y.; Liu, X.P.; Li, K.; Li, X. A maximum entropy method to extract urban land by combining MODIS reflectance, MODIS NDVI, and DMSP-OLS data. Int. J. Remote Sens. 2014, 35, 6708–6727. [Google Scholar] [CrossRef]
  25. Hu, Y.; Li, H.; Wu, D.; Chen, W.; Zhao, X.; Hou, M.; Li, A.; Zhu, Y. LAI-indicated vegetation dynamic in ecologically fragile region: A case study in the Three-North Shelter Forest program region of China. Ecol. Indic. 2021, 120, 106932. [Google Scholar] [CrossRef]
  26. Zheng, K.; Tan, L.; Sun, Y.; Wu, Y.; Duan, Z.; Xu, Y.; Gao, C. Impacts of climate change and anthropogenic activities on vegetation change: Evidence from typical areas in China. Ecol. Indic. 2021, 126, 107648. [Google Scholar] [CrossRef]
  27. Zhu, Z.; Piao, S.; Myneni, R.B.; Huang, M.; Zeng, Z.; Canadell, J.G.; Ciais, P.; Sitch, S.; Friedlingstein, P.; Arneth, A.; et al. Greening of the Earth and its drivers. Nat. Clim. Chang. 2016, 6, 791–795. [Google Scholar] [CrossRef]
  28. Li, H.X.; Liu, G.H.; Fu, B.J. Spatial variations of rain-use efficiency along a climate gradient on the Tibetan Plateau: A satellite-based analysis. Int. J. Remote Sens. 2013, 34, 7487–7503. [Google Scholar] [CrossRef]
  29. Senanayake, S.; Pradhan, B.; Huete, A.; Brennan, J. Spatial modeling of soil erosion hazards and crop diversity change with rainfall variation in the Central Highlands of Sri Lanka. Sci. Total Environ. 2022, 806, 150405. [Google Scholar] [CrossRef]
  30. Li, C.; Zhao, J.; Shi, Y.; Hu, X. The impact of human activities on net primary productivity based on the coefficient of variation: A case study of the Shiyang River Basin. Acta Ecol. Sin. 2016, 36, 4034–4044. [Google Scholar]
  31. Qi, X.; Jia, J.; Liu, H.; Lin, Z. Relative importance of climate change and human activities for vegetation changes on China’s silk road economic belt over multiple timescales. Catena 2019, 180, 224–237. [Google Scholar] [CrossRef]
  32. Kastner, T.; Matej, S.; Forrest, M.; Gingrich, S.; Haberl, H.; Hickler, T.; Krausmann, F.; Lasslop, G.; Niedertscheider, M.; Plutzar, C.; et al. Land use intensification increasingly drives the spatiotemporal patterns of the global human appropriation of net primary production in the last century. Glob. Chang. Biol. 2022, 28, 307–322. [Google Scholar] [CrossRef] [PubMed]
  33. Zhou, X.; Yamaguchi, Y.; Arjasakusuma, S. Distinguishing the vegetation dynamics induced by anthropogenic factors using vegetation optical depth and AVHRR NDVI: A cross-border study on the Mongolian Plateau. Sci. Total Environ. 2018, 616–617, 730–743. [Google Scholar] [CrossRef] [PubMed]
  34. Qu, S.; Wang, L.; Lin, A.; Zhu, H.; Yuan, M. What drives the vegetation restoration in Yangtze River basin, China: Climate change or anthropogenic factors? Ecol. Indic. 2018, 90, 438–450. [Google Scholar] [CrossRef]
  35. Wu, J.; Li, M.; Zhang, X.; Fiedler, S.; Gao, Q.; Zhou, Y.; Cao, W.; Hassan, W.; Margarint, M.C.; Tarolli, P.; et al. Disentangling climatic and anthropogenic contributions to nonlinear dynamics of alpine grassland productivity on the Qinghai-Tibetan Plateau. J. Environ. Manag. 2021, 281, 111875. [Google Scholar] [CrossRef]
  36. Yuan, W.P.; Zheng, Y.; Piao, S.L.; Ciais, P.; Lombardozzi, D.; Wang, Y.P.; Ryu, Y.; Chen, G.X.; Dong, W.J.; Hu, Z.M.; et al. Increased atmospheric vapor pressure deficit reduces global vegetation growth. Sci. Adv. 2019, 5, eaax1396. [Google Scholar] [CrossRef] [Green Version]
  37. Zhang, Y.; Zhang, C.; Wang, Z.; Chen, Y.; Gang, C.; An, R.; Li, J. Vegetation dynamics and its driving forces from climate change and human activities in the Three-River Source Region, China from 1982 to 2012. Sci. Total Environ. 2016, 563–564, 210–220. [Google Scholar] [CrossRef]
  38. Dong, Q.-M.; Zhao, X.-Q.; Wu, G.-L.; Shi, J.-J.; Ren, G.-H. A review of formation mechanism and restoration measures of “black-soil-type” degraded grassland in the Qinghai-Tibetan Plateau. Environ. Earth Sci. 2013, 70, 2359–2370. [Google Scholar] [CrossRef]
  39. Piao, S.; Ciais, P.; Friedlingstein, P.; Peylin, P.; Reichstein, M.; Luyssaert, S.; Margolis, H.; Fang, J.; Barr, A.; Chen, A.; et al. Net carbon dioxide losses of northern ecosystems in response to autumn warming. Nature 2008, 451, 49–52. [Google Scholar] [CrossRef]
  40. Bai, Y.; Guo, C.; Degen, A.A.; Ahmad, A.A.; Wang, W.; Zhang, T.; Li, W.; Ma, L.; Huang, M.; Zeng, H.; et al. Climate warming benefits alpine vegetation growth in Three-River Headwater Region, China. Sci. Total Environ. 2020, 742, 140574. [Google Scholar] [CrossRef]
  41. Fu, Z.; Ciais, P.; Prentice, I.C.; Gentine, P.; Makowski, D.; Bastos, A.; Luo, X.; Green, J.K.; Stoy, P.C.; Yang, H.; et al. Atmospheric dryness reduces photosynthesis along a large range of soil water deficits. Nat. Commun. 2022, 13, 989. [Google Scholar] [CrossRef]
  42. Xiu, L.; Yao, X.; Chen, M.; Yan, C. Effect of Ecological Construction Engineering on Vegetation Restoration: A Case Study of the Loess Plateau. Remote Sens. 2021, 13, 1407. [Google Scholar] [CrossRef]
  43. Shao, Q.; Cao, W.; Fan, J.; Huang, L.; Xu, X. Effects of an ecological conservation and restoration project in the Three-River Source Region, China. J. Geogr. Sci. 2016, 27, 183–204. [Google Scholar] [CrossRef]
  44. Shen, X.; An, R.; Feng, L.; Ye, N.; Zhu, L.; Li, M. Vegetation changes in the Three-River Headwaters Region of the Tibetan Plateau of China. Ecol. Indic. 2018, 93, 804–812. [Google Scholar] [CrossRef]
  45. Liu, Y.; Liu, R.; Chen, J.M. Retrospective retrieval of long-term consistent global leaf area index (1981–2011) from combined AVHRR and MODIS data. J. Geophys. Res.-Biogeosci. 2012, 117, 1–14. [Google Scholar] [CrossRef]
  46. Deng, F.; Chen, J.M.; Plummer, S.; Chen, M.Z.; Pisek, J. Algorithm for global leaf area index retrieval using satellite imagery. IEEE Trans. Geosci. Remote Sens. 2006, 44, 2219–2229. [Google Scholar] [CrossRef] [Green Version]
  47. Tucker, C.J.; Pinzon, J.E.; Brown, M.E.; Slayback, D.A.; Pak, E.W.; Mahoney, R.; Vermote, E.F.; El Saleous, N. An extended AVHRR 8-km NDVI dataset compatible with MODIS and SPOT vegetation NDVI data. Int. J. Remote Sens. 2005, 26, 4485–4498. [Google Scholar] [CrossRef]
  48. Li, X.; Lu, H.; Yu, L.; Yang, K. Comparison of the Spatial Characteristics of Four Remotely Sensed Leaf Area Index Products over China: Direct Validation and Relative Uncertainties. Remote Sens. 2018, 10, 148. [Google Scholar] [CrossRef] [Green Version]
  49. Zhang, X.R.; Zhu, L.; Sun, H.F.; Chu, S.S. Validation and inter-comparison of the FY-3B/MERSI LAI product with GLOBMAP and MYD15A2H. Int. J. Remote Sens. 2020, 41, 9256–9282. [Google Scholar] [CrossRef]
  50. Jiang, C.; Ryu, Y.; Fang, H.; Myneni, R.; Claverie, M.; Zhu, Z. Inconsistencies of interannual variability and trends in long-term satellite leaf area index products. Glob. Chang. Biol. 2017, 23, 4133–4146. [Google Scholar] [CrossRef]
  51. Peng, S.; Ding, Y.; Liu, W.; Li, Z. 1 km monthly temperature and precipitation dataset for China from 1901 to 2017. Earth Syst. Sci. Data 2019, 11, 1931–1946. [Google Scholar] [CrossRef] [Green Version]
  52. Abatzoglou, J.T.; Dobrowski, S.Z.; Parks, S.A.; Hegewisch, K.C. TerraClimate, a high-resolution global dataset of monthly climate and climatic water balance from 1958–2015. Sci. Data 2018, 5, 170191. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Zhai, X.; Liang, X.; Yan, C.; Xing, X.; Jia, H.; Wei, X.; Feng, K. Vegetation Dynamic Changes and Their Response to Ecological Engineering in the Sanjiangyuan Region of China. Remote Sens. 2020, 12, 4035. [Google Scholar] [CrossRef]
  54. Li, Y.; Zheng, Z.; Qin, Y.; Rong, P. Relative contributions of natural and man-made factors to vegetation cover change of environmentally sensitive and vulnerable areas of China. J. Clean Prod. 2021, 321, 128917. [Google Scholar] [CrossRef]
  55. Evans, J.; Geerken, R. Discrimination between climate and human-induced dryland degradation. J. Arid. Environ. 2004, 57, 535–554. [Google Scholar] [CrossRef]
  56. Wessels, K.J.; Prince, S.D.; Malherbe, J.; Small, J.; Frost, P.E.; VanZyl, D. Can human-induced land degradation be distinguished from the effects of rainfall variability? A case study in South Africa. J. Arid. Environ. 2007, 68, 271–297. [Google Scholar] [CrossRef]
  57. Liu, X.; Zhang, J.; Zhu, X.; Pan, Y.; Liu, Y.; Zhang, D.; Lin, Z. Spatiotemporal changes in vegetation coverage and its driving factors in the Three-River Headwaters Region during 2000–2011. J. Geogr. Sci. 2014, 24, 288–302. [Google Scholar] [CrossRef]
  58. Deng, L.S.; Liu, G.Z.; Wu, X.; Wang, Y.P.; Tong, M.H.; Liu, B.; Wang, K.D.; Peng, Y.T.; Kong, X.P. Adipose Derived Mesenchymal Stem Cells Efficiently Rescue Carbon Tetrachloride-Induced Acute Liver Failure in Mouse. Sci. World J. 2014, 2014, 103643. [Google Scholar] [CrossRef] [Green Version]
  59. Chen, Y.P.; Wang, K.B.; Lin, Y.S.; Shi, W.Y.; Song, Y.; He, X.H. Balancing green and grain trade. Nat. Geosci. 2015, 8, 739–741. [Google Scholar] [CrossRef]
  60. Han, Z.M.; Huang, S.Z.; Huang, Q.; Bai, Q.J.; Leng, G.Y.; Wang, H.; Zhao, J.; Wei, X.T.; Zheng, X.D. Effects of vegetation restoration on groundwater drought in the Loess Plateau, China. J. Hydrol. 2020, 591, 125566. [Google Scholar] [CrossRef]
  61. Dong, S.; Shang, Z.; Gao, J.; Boone, R.B. Enhancing sustainability of grassland ecosystems through ecological restoration and grazing management in an era of climate change on Qinghai-Tibetan Plateau. Agric. Ecosyst. Environ. 2020, 287, 106684. [Google Scholar] [CrossRef]
  62. Shi, Y.; Jin, N.; Ma, X.; Wu, B.; He, Q.; Yue, C.; Yu, Q. Attribution of climate and human activities to vegetation change in China using machine learning techniques. Agric. For. Meteorol. 2020, 294, 108146. [Google Scholar] [CrossRef]
  63. He, P.; Ma, J.; Han, Z.; Shi, M.; Xu, D.; Sun, Z. Relationship between multi-scale climate factors and performance of ecological engineering on the Loess Plateau, China. J. For. Res. 2021, 33, 789–800. [Google Scholar] [CrossRef]
  64. Corney, P.M.; Le Duc, M.G.; Smart, S.M.; Kirby, K.J.; Bunce, R.G.H.; Marrs, R.H. The effect of landscape-scale environmental drivers on the vegetation composition of British woodlands. Biol. Conserv. 2004, 120, 491–505. [Google Scholar] [CrossRef]
  65. Bao, Z.X.; Zhang, J.Y.; Wang, G.Q.; Guan, T.S.; Jin, J.L.; Liu, Y.L.; Li, M.; Ma, T. The sensitivity of vegetation cover to climate change in multiple climatic zones using machine learning algorithms. Ecol. Indic. 2021, 124, 107443. [Google Scholar] [CrossRef]
  66. Guo, B.; Wang, J.; Mantravadi, V.S.; Zhang, L.; Liu, G. Effect of climate and ecological restoration on vegetation changes in the “Three-River Headwaters” region based on remote sensing technology. Environ. Sci. Pollut. Res. 2021, 29, 16436–16448. [Google Scholar] [CrossRef]
  67. Cong, N.; Shen, M.G.; Zu, J.X.; Zhang, Y.J. Spatial sampling inconsistency leads to differences in phenological sensitivity to warming between natural and experiment sites. Sci. Bull. 2019, 64, 961–963. [Google Scholar] [CrossRef] [Green Version]
  68. Berdimbetov, T.; Ma, Z.-G.; Shelton, S.; Ilyas, S.; Nietullaeva, S. Identifying Land Degradation and its Driving Factors in the Aral Sea Basin From 1982 to 2015. Front. Earth Sci. 2021, 9, 690000. [Google Scholar] [CrossRef]
  69. Kong, D.; Zhang, Q.; Singh, V.P.; Shi, P. Seasonal vegetation response to climate change in the Northern Hemisphere (1982–2013). Glob. Planet. Chang. 2017, 148, 1–8. [Google Scholar] [CrossRef] [Green Version]
  70. Fang, J.; Piao, S.; He, J.; Ma, W. Increasing terrestrial vegetation activity in China, 1982–1999. Sci. China Ser. C-Life Sci. 2004, 47, 229–240. [Google Scholar]
  71. Chuai, X.W.; Huang, X.J.; Wang, W.J.; Bao, G. NDVI, temperature and precipitation changes and their relationships with different vegetation types during 1998–2007 in Inner Mongolia, China. Int. J. Climatol. 2013, 33, 1696–1706. [Google Scholar] [CrossRef]
  72. Pan, X.; Song, Y.B.; Liu, G.F.; Hu, Y.K.; Ye, X.H.; Cornwell, W.K.; Prinzing, A.; Dong, M.; Cornelissen, J.H.C. Functional traits drive the contribution of solar radiation to leaf litter decomposition among multiple arid-zone species. Sci. Rep. 2015, 5, 13217. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Qu, S.; Wang, L.; Lin, A.; Yu, D.; Yuan, M. Distinguishing the impacts of climate change and anthropogenic factors on vegetation dynamics in the Yangtze River Basin, China. Ecol. Indic. 2020, 108, 105724. [Google Scholar] [CrossRef]
  74. Zhang, X.; Jin, X. Vegetation dynamics and responses to climate change and anthropogenic activities in the Three-River Headwaters Region, China. Ecol. Indic. 2021, 131, 108223. [Google Scholar] [CrossRef]
  75. Zhang, X.; Yue, Y.; Tong, X.; Wang, K.; Qi, X.; Deng, C.; Brandt, M. Eco-engineering controls vegetation trends in southwest China karst. Sci. Total Environ. 2021, 770, 145160. [Google Scholar] [CrossRef]
  76. Li, S.; Yang, S.; Liu, X.; Liu, Y.; Shi, M. NDVI-Based Analysis on the Influence of Climate Change and Human Activities on Vegetation Restoration in the Shaanxi-Gansu-Ningxia Region, Central China. Remote Sens. 2015, 7, 11163–11182. [Google Scholar] [CrossRef] [Green Version]
  77. Liu, L.; Shao, Q.; Cao, W.; Wu, D.; Gong, G.; Fan, J. Cost-Benefit Analysis of the Ecological Protects in the Three-River Headwaters Region Based on Ecosystem Services Values. Acta Agrestia Sin. 2018, 26, 30–39. [Google Scholar]
  78. Fang, Y. Managing the three-rivers headwater region, china: From ecological engineering to social engineering. Ambio 2013, 42, 566–576. [Google Scholar] [CrossRef] [Green Version]
  79. Fan, J.W.; Shao, Q.Q.; Liu, J.Y.; Wang, J.B.; Harris, W.; Chen, Z.Q.; Zhong, H.P.; Xu, X.L.; Liu, R.G. Assessment of effects of climate change and grazing activity on grassland yield in the Three Rivers Headwaters Region of Qinghai-Tibet Plateau, China. Environ. Monit. Assess. 2010, 170, 571–584. [Google Scholar] [CrossRef]
  80. Tong, X.; Brandt, M.; Yue, Y.; Horion, S.; Wang, K.; Keersmaecker, W.D.; Tian, F.; Schurgers, G.; Xiao, X.; Luo, Y.; et al. Increased vegetation growth and carbon stock in China karst via ecological engineering. Nat. Sustain. 2018, 1, 44–50. [Google Scholar] [CrossRef]
  81. Xie, S.; Mo, X.; Hu, S.; Liu, S. Contributions of climate change, elevated atmospheric CO2 and human activities to ET and GPP trends in the Three-North Region of China. Agric. For. Meteorol. 2020, 295, 108183. [Google Scholar] [CrossRef]
  82. Tang, L.; Chen, X.; Cai, X.; Li, J. Disentangling the roles of land-use-related drivers on vegetation greenness across China. Environ. Res. Lett. 2021, 16, 124033. [Google Scholar] [CrossRef]
Figure 1. Location and topographical features (a) and land cover types (b) of the study area.
Figure 1. Location and topographical features (a) and land cover types (b) of the study area.
Remotesensing 14 02386 g001
Figure 2. The trends in the LAI values over the LP and TRSR for the entire period (1982–2019); two stages were set according to the initiation timing of the ecological engineering project (1982–1999 and 2000–2019 for the LP; 1982–2004 and 2005–2019 for the TRSR).
Figure 2. The trends in the LAI values over the LP and TRSR for the entire period (1982–2019); two stages were set according to the initiation timing of the ecological engineering project (1982–1999 and 2000–2019 for the LP; 1982–2004 and 2005–2019 for the TRSR).
Remotesensing 14 02386 g002
Figure 3. Spatial distributions of the LAI trends over the study areas during the three periods: 1982–2019 (a), before the implementation of the ecological engineering project (b), and after the implementation of the ecological engineering project (c). The inset panels show the pixels where LAI trends are statistically significant at p < 0.05.
Figure 3. Spatial distributions of the LAI trends over the study areas during the three periods: 1982–2019 (a), before the implementation of the ecological engineering project (b), and after the implementation of the ecological engineering project (c). The inset panels show the pixels where LAI trends are statistically significant at p < 0.05.
Remotesensing 14 02386 g003
Figure 4. Interannual changes in annual mean climatic variables (TMP, PRE, RAD, and VPD) during 1982–2019. TMP, PRE, RAD, and VPD represent temperature, precipitation, downward surface shortwave radiation, and vapor pressure deficit, respectively.
Figure 4. Interannual changes in annual mean climatic variables (TMP, PRE, RAD, and VPD) during 1982–2019. TMP, PRE, RAD, and VPD represent temperature, precipitation, downward surface shortwave radiation, and vapor pressure deficit, respectively.
Remotesensing 14 02386 g004
Figure 5. Spatial patterns of the trends in annual mean climatic variables (TMP, PRE, RAD, and VPD) for the three periods: 1982–2019 (a1d1), before the implementation of ecological engineering project (a2d2), after the implementation of ecological engineering project (a3d3). TMP, PRE, RAD, and VPD represent temperature, precipitation, downward surface shortwave radiation, and vapor pressure deficit, respectively.
Figure 5. Spatial patterns of the trends in annual mean climatic variables (TMP, PRE, RAD, and VPD) for the three periods: 1982–2019 (a1d1), before the implementation of ecological engineering project (a2d2), after the implementation of ecological engineering project (a3d3). TMP, PRE, RAD, and VPD represent temperature, precipitation, downward surface shortwave radiation, and vapor pressure deficit, respectively.
Remotesensing 14 02386 g005
Figure 6. Spatial patterns of the correlation coefficients between LAI and TMP (first column), PRE (second column), RAD (third column), VPD (fourth column), and spatial patterns of dominant climatic constraint on LAI (fifth column). They were estimated for 1982–2019 (a1e1), before the implementation of the ecological engineering project (a2e2), after the implementation of ecological engineering (a3e3), respectively. TMP, PRE, RAD, and VPD represent temperature, precipitation, downward surface shortwave radiation, and vapor pressure deficit, respectively.
Figure 6. Spatial patterns of the correlation coefficients between LAI and TMP (first column), PRE (second column), RAD (third column), VPD (fourth column), and spatial patterns of dominant climatic constraint on LAI (fifth column). They were estimated for 1982–2019 (a1e1), before the implementation of the ecological engineering project (a2e2), after the implementation of ecological engineering (a3e3), respectively. TMP, PRE, RAD, and VPD represent temperature, precipitation, downward surface shortwave radiation, and vapor pressure deficit, respectively.
Remotesensing 14 02386 g006
Figure 7. Spatial patterns of the linear regression coefficients between the LAI and TMP (a), PRE (b), RAD (c), VPD (d) during 1982–2019.
Figure 7. Spatial patterns of the linear regression coefficients between the LAI and TMP (a), PRE (b), RAD (c), VPD (d) during 1982–2019.
Remotesensing 14 02386 g007
Figure 8. Spatial patterns of the relative contribution of human activities to vegetation greening (p < 0.05) during the three periods: 1982–2019 (a), before the implementation of ecological engineering (b), and after the implementation of ecological engineering (c).
Figure 8. Spatial patterns of the relative contribution of human activities to vegetation greening (p < 0.05) during the three periods: 1982–2019 (a), before the implementation of ecological engineering (b), and after the implementation of ecological engineering (c).
Remotesensing 14 02386 g008
Table 1. Determination criteria of driving factors on vegetation greening and calculation methods of their contribution.
Table 1. Determination criteria of driving factors on vegetation greening and calculation methods of their contribution.
Vegetation GreeningDetermination CriteriaContribution (%)
L A I a c t u a l   Trend L A I p r e   Trend L A I r e s   Trend Climate VariabilitiesHuman Activities
>0>0>0 S l o p e ( L A I p r e ) S l o p e ( L A I a c t u a l ) S l o p e ( L A I r e s ) S l o p e ( L A I a c t u a l )
<0>00100
>0<01000
<0<0<0 S l o p e ( L A I p r e ) S l o p e ( L A I a c t u a l ) S l o p e ( L A I r e s ) S l o p e ( L A I a c t u a l )
<0>01000
>0<00100
Table 2. Summary of LAI trends in the study regions.
Table 2. Summary of LAI trends in the study regions.
RegionsPeriodsIncreased Area (%)Significantly Increased Area (%)Decreased Area (%)Significantly Decreased Area (%)
LP1982–201997.591.42.50.5
1982–199980.323.319.70.9
2000–201998.994.71.10.3
TRSR1982–201986.065.314.13.5
1982–200466.226.433.83.9
2005–201970.517.029.52.6
Table 3. Summary of the correlations between LAI and TMP, PRE, RAD, and VPD in the study regions.
Table 3. Summary of the correlations between LAI and TMP, PRE, RAD, and VPD in the study regions.
BetweenRelationshipLPTRSR
1982–20191982–19992000–20191982–20191982–20042005–2019
LAI-TMPPositive area (%)96.6571.5288.6290.0581.2570.61
Significant positive area (%)72.098.996.5057.0917.563.55
Negative area (%)3.3528.4811.389.9518.7529.39
Significant negative area (%)0.220.320.010.800.250.41
LAI-PREPositive area (%)95.2471.9595.0088.1141.8282.57
Significant positive area (%)42.376.6735.1730.941.2515.46
Negative area (%)4.7628.055.0011.8958.1817.43
Significant negative area (%)0.191.350.040.251.810.17
LAI-RADPositive area (%)55.2636.5256.9811.7145.7215.54
Significant positive area (%)9.671.802.700.080.980.17
Negative area (%)44.7463.4843.0288.2954.2884.46
Significant negative area (%)11.023.820.3620.642.0312.58
LAI-VPDPositive area (%)80.2447.4316.6626.7121.9045.50
Significant positive area (%)26.984.010.190.310.462.36
Negative area (%)19.7652.5783.3473.2978.1054.50
Significant negative area (%)1.694.1312.0113.4410.511.33
Table 4. Transfer matrix of different land cover types from 1980 to 2018 in the study regions (unit: km2).
Table 4. Transfer matrix of different land cover types from 1980 to 2018 in the study regions (unit: km2).
PeriodTypeCroplandForestShrublandGrasslandWaterbodiesUrbanDesert
1980–2000
LP
Cropland201,270.44148.8093.413364.03792.368.89608.26
Forest208.4952,208.5660.11322.9131.450.0098.74
Shrubland40.81122.6438,637.24344.418.150.0066.18
Grassland1077.27385.56424.28254,937.44389.212.622103.13
Waterbodies255.486.1610.59287.907948.660.00129.21
Urban1235.3237.495.01242.1334.6413,126.3232.19
Desert319.478.1843.281568.61243.710.0040,431.11
2000–2018
LP
Cropland117,484.226639.214311.6154,293.922315.977088.912599.08
Forest7993.5930,779.403830.0211,692.15238.07285.55566.68
Shrubland5051.893719.1919,705.7510,686.85119.6392.60325.40
Grassland57,769.819927.4410,364.27164,034.371857.111972.0912,574.92
Waterbodies2391.27263.56154.342028.062921.85231.08652.60
Urban13,686.50918.06383.515605.51539.484870.88919.56
Desert1604.24403.07287.9910,301.00591.73151.9524,717.00
1980–2005
TRSR
Cropland745.284.250.0086.002.340.002.87
Forest0.003128.34.1421.831.000.001.03
Shrubland0.004.7711,037.21.271.970.000.00
Grassland3.2418.8037.51237,798.15312.442.33192.60
Waterbodies2.560.002.6268.1416,626.610.0094.05
Urban4.350.000.004.661.8261.800.00
Desert2.590.002.38924.35485.080.0078,932.47
2005–2018
TRSR
Cropland456.3119.0029.67335.2129.278.0629.52
Forest23.071578.30194.961257.7123.960.4094.47
Shrubland18.58194.494763.105716.5349.713.25221.99
Grassland284.391318.765828.42210,490.874350.7336.9330,601.93
Waterbodies32.6023.1344.434524.1010,457.732.782470.63
Urban17.603.126.4768.644.9916.976.22
Desert2.7720.32159.9515,382.961793.671.6246,636.25
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gao, J.; Zhang, Y.; Zheng, Z.; Cong, N.; Zhao, G.; Zhu, Y.; Chen, Y.; Sun, Y.; Zhang, J.; Zhang, Y. Ecological Engineering Projects Shifted the Dominance of Human Activity and Climate Variability on Vegetation Dynamics. Remote Sens. 2022, 14, 2386. https://doi.org/10.3390/rs14102386

AMA Style

Gao J, Zhang Y, Zheng Z, Cong N, Zhao G, Zhu Y, Chen Y, Sun Y, Zhang J, Zhang Y. Ecological Engineering Projects Shifted the Dominance of Human Activity and Climate Variability on Vegetation Dynamics. Remote Sensing. 2022; 14(10):2386. https://doi.org/10.3390/rs14102386

Chicago/Turabian Style

Gao, Jie, Yangjian Zhang, Zhoutao Zheng, Nan Cong, Guang Zhao, Yixuan Zhu, Yao Chen, Yihan Sun, Jianshuang Zhang, and Yu Zhang. 2022. "Ecological Engineering Projects Shifted the Dominance of Human Activity and Climate Variability on Vegetation Dynamics" Remote Sensing 14, no. 10: 2386. https://doi.org/10.3390/rs14102386

APA Style

Gao, J., Zhang, Y., Zheng, Z., Cong, N., Zhao, G., Zhu, Y., Chen, Y., Sun, Y., Zhang, J., & Zhang, Y. (2022). Ecological Engineering Projects Shifted the Dominance of Human Activity and Climate Variability on Vegetation Dynamics. Remote Sensing, 14(10), 2386. https://doi.org/10.3390/rs14102386

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop